clear variables
close all
Racunanje kompleksne struje generatora i snage na ulazu voda
km = 1e3;
Rp = 0.08/km;
Lp = 1.34e-3/km;
Cp = 8.6e-9/km;
Gp = 37.5e-9/km;
D = 900*km;
f = 50;
P = 100e6;
k = 1;
U = 220e3;
Zc = sqrt((Rp + 1i*2*pi*f*Lp) / (Gp + 1i*2*pi*f*Cp))
Zc = 3.9673e+02 - 3.4586e+01i
gamma = sqrt((Rp + 1i*2*pi*f*Lp) * (Gp + 1i*2*pi*f*Cp))
gamma = 1.0832e-07 + 1.0706e-06i
a11 = cosh(gamma*D)
a11 = 0.5734 + 0.0802i
a12 = Zc*sinh(gamma*D)
a12 = 5.0642e+01 + 3.2541e+02i
a21 = (1/Zc)*sinh(gamma*D)
a21 = -0.0000 + 0.0021i
a22 = cosh(gamma*D)
a22 = 0.5734 + 0.0802i
ID = P/(U*k)
ID = 454.5455
U0 = a11*U + a12*ID
U0 = 1.4916e+05 + 1.6556e+05i
I0 = a21*U + a22*ID
I0 = 2.5169e+02 + 4.9322e+02i
S0 = U0*conj(I0)
S0 = 1.1920e+08 - 3.1899e+07i
P0 = real(S0)
P0 = 1.1920e+08
Odredjivanje zamenske Pi sheme
syms Zc gamma D Z1 Z2 a11 a12 a21 a22 U0 I0 UD ID
assume(Z1~=0 & Z2~=0 & a11~=0 & a12~=0)
jednacine = [U0 == Z2*(ID+UD/Z1) + UD,...
I0 == U0/Z1 + (U0-UD)/Z2]
jednacine = 
resenje = solve(jednacine, [U0,I0])
resenje = struct with fields:
U0: (UD*Z1 + UD*Z2 + ID*Z1*Z2)/Z1 I0: (2*UD*Z1 + UD*Z2 + ID*Z1^2 + ID*Z1*Z2)/Z1^2
jed1 = a11 == simplify(subs(resenje.U0/UD, ID,0))
jed1 = 
jed2 = a12 == subs(resenje.U0/ID, UD,0)
jed2 = 
resenjeZ = solve([jed1, jed2], [Z1,Z2])
resenjeZ = struct with fields:
Z1: a12/(a11 - 1) Z2: a12
Z1 = subs(resenjeZ.Z1, [a11 a12], [cosh(gamma*D) Zc*sinh(gamma*D)])
Z1 = 
Z2 = subs(resenjeZ.Z2, [a11 a12], [cosh(gamma*D) Zc*sinh(gamma*D)])
Z2 =